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We consider a spin-1/2 anisotropic XY model with time-dependent spin-spin coupling as means of 
creating long-distance entanglement. We predict the emergence of significant entanglement between 
the first and the last spin whenever the ac part of the coupling has a frequency matching the Zee- 
man splitting. In particular, we find that the concurrence assumes its maximum with a vanishing 
dc part. Mapping the time-dependent Hamiltonian within a rotating-wave approximation to an 
effective static model provides qualitative and quantitative understanding of this entanglement res- 
onance. Numerical results for the duration of the entanglement creation and its length dependence 
substantiate the effective static picture. 
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I. INTRODUCTION 

Entanglement is a key resource for many quantum in- 
formation and computation protocols, such as teleporta- 
tion jT] , superdense coding [5] , and cryptography [3] . The 
successful storage and transfer of quantum information 
requires effective mechanisms to create entangled states 
over large distances. Since entanglement is generated 
mostly by local interactions, it is initially short-ranged 
and, thus, has to be distributed via quantum channels. 
Lately it has been noticed that spin chains are promising 
candidates for this task Various spin-spin interac- 
tions, like e.g. Ising or Heisenberg coupling, have been 
considered for entanglement creation, and their static as 
well as their dynamical properties have been investigated 
[3 El [II [3 [9] . Spin chains thus turned out to be efficient 
quantum channels for controlled entanglement distribu- 
tion. 

A particular spin chain is the quantum anisotropic XY 
model. Irrespective of the magnitude of the anisotropy, 
it can be solved exactly with the help of a Jordan- 
Wigner transformation and therefore became a paradig- 
matic model in many-body physics jlOj . In the context 
of quantum information [TT] , its experimental implemen- 
tation with optical lattices [T2], quantum dots [T3], and 
Josephson junctions [14j has been proposed. 

Thus far, most studies consider transfer of entangle- 
ment rather than its generation |15j . or its presence in 
systems with static interactions [HI [171 [E] • In this pa- 
per, by contrast, we analyze a spin-1/2 XY chain with pe- 
riodically time- dependent nearest-neighbor coupling with 
separable initial state and find entanglement creation 
between the first and the last spin of the chain. As 
it may be difficult to access individual spins in a con- 
trolled manner, we restrict ourselves to chains with global 
time-dependent spin-spin coupling. Remarkably, entan- 
glement created in that way turns out to be significantly 
larger than the one in related static systems [THl [20] . We 
gain further insight by mapping the time-dependent spin 
chain to a static model, which also provides information 
on the length dependence and the duration of the entan- 
glement creation. 



II. SPIN CHAIN HAMILTONIAN 

The anisotropic XY model in a transverse field B and 
with time-dependent nearest-neighbor coupling J(t) is 
described by the Hamiltonian form (we put h — 1) 

^= yE<+^E [(1+7)<^,1+i+(1-7)«+i], 

11—1 n—1 

(1) 

where the cr's are the usual Pauli matrices and 7 de- 
notes the anisotropy parameter. We focus on situations 
in which the coupling strength J(t) is smaller than the 
field strength B and are interested in the spectral re- 
sponse of the chain when the coupling is periodically 
modulated. All other parameters have arbitrary but fixed 
values. We also suppose that the spins are initially un- 
coupled, J(t) = for < < 0, and cooled down to the 
fully-aligned separable state 

IV'(i-O)) = 10000 ...) , (2) 

which is the ground state of the Hamiltonian ([l]) with 
J = 0. At t = 0, we switch on a coupling consisting 
of a dc contribution Jq and a sinusoidal ac part with 
amplitude Ji, 

J{t > 0) = Jo + Ji sin(tJdi) . (3) 

In the limit Ji — > 0, the coupling suddenly switches to a 
constant value, while for Ji ^ 0, we are able to probe the 
frequency-dependent response of the system. We quan- 
tify entanglement between the two ends of the chain with 
the help of the concurrence C = max{Ai — A2 — A3 — A4, 0}. 
The A's are the ordered square roots of the eigenvalues of 
p{ay ® ^y)P*{'^y ® ^y) ^^th p being the reduced density 
matrix of the two spins |21j . 

III. ENTANGLEMENT RESONANCE 

By direct numerical integration, we investigated the 
time evolution of the concurrence between the first 
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FIG. 1: (color online) Maximum concurrence obtained in a time window up to 4A''/max(Jo, Ji), between spins 1 and A'^ for 
different frequencies and chain lengths A'' with Ji — O.IB, 7 = 1, and Jo = (solid), O.OIB (dashed), O.IB (dash-dotted). The 
effective anisotropy thus has the values 7 = 7J1/2J0 — 00, 5, 0.5. 



and the last spin for different driving frequencies lu^, 
chain lengths N, and parameters 7, Ji, Jq. We de- 
termined the maximal concurrence in the time interval 
[0, . . . ,4A^/max(Jo, Ji)]. The results shown in Fig. [l]re- 
veal that at cod = 2B, irrespective of the other parame- 
ters, the concurrence assumes during that time interval 
a value close to unity and is significantly larger than for 
other frequencies — we term this entanglement resonance. 
Height and width of the resonance peak depend on the 
intensities Jq and Ji and on the chain length N; see 
Fig. |2] We also notice the existence of a much smaller 
secondary peak at uj^ = B. Contrary to the main peak, 
its amplitude strongly decreases with decreasing coupling 
intensity and increasing chain length. We henceforth fo- 
cus on the peak at tOd — 2B. 



A. Rotating-wave approximation (RWA) 

Deeper understanding of the observed entanglement 
resonance can be gained by analyzing the time-dependent 
Hamiltonian ([T]) within rotating-wave approximation. 
Since entanglement properties of a system are not 
changed by local unitary operations on individual subsys- 
tems, it is convenient to transform the XY Hamiltonian 
to the interaction picture, H = exp(iHQt)H cxp(~iHQt), 
with Hq — (-B/2) '^i- By introducing the shift opera- 
tors cr'' 



\{<y^ ± ic^), we obtain 
J{t) 



N 



(4) 



2iBt 



+7e-"V-<+i+7e 



n '-'n+1 



The first two terms swap the states of spins n and n -I- 1, 
while the last two terms pairwise create (destroy) excita- 
tions, which here is the origin of entanglement generation. 
If the driving frequency obeys the resonance condition 
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FIG. 2: (color online) Maximum obtained concurrence be- 
tween the end spins for ujd — 2B as a function of the dc inter- 
action Jo for various chain lengths, Ji = O.IB and 7 = 1. The 
solid lines are obtained with the full time-dependent Hamil- 
tonian while the dashed lines mark the RWA solution. 



LJd — 2B and, moreover, is much larger than both Jq and 
Ji, we can within RWA replace the Hamiltonian (|4| by 
its time average 



^ r 



n=l 



(5) 



with the effective anisotropy 7 = 7J1/2J0. This means 
that for resonant driving, the time dependent XY model 
([ij can be mapped to the static XY model ([5| without 
any Zeeman field. In both cases, the entanglement gener- 
ated between the two end spins is maximal and controlled 
by the parameter 7 and the chain length N. Note that 
Jo — > corresponds to the infinitely anisotropic limit 
7 — > 00. Figure [2] shows the concurrence between the 
end spins as a function of the anisotropy parameter 7 for 
resonant driving. Two important points are worth be- 
ing mentioned: First, the concurrence approaches unity 
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FIG. 3; (color online) Entanglement dynamics for a chain of 
8 spins with driving frequency cud ~ 2B for two different ef- 
fective anisotropies 7 = 7J1/2J0. As in Fig. [1] Ji = O.IB, 
7 = 1. Solid lines mark the exact numerical solution, while 
the dashed lines are computed within RWA. The symbols 
mark the time evolution for switching off the driving (squares) 
and for changing the frequency to u = 3B (circles) after the 
concurrence maximum is reached at time tarrivai . Both results 
cannot be distinguished for the chosen resolution. 

in the limit of vanishing Jq, i.e. for infinite 7. In this 
limit, the amount of entanglement no longer depends on 
Ji and 7. Second, the agreement of the exactly eval- 
uated concurrence and the RWA solution is excellent, 
which demonstrates that RWA is appropriate. Moreover, 
Fig. [3] shows that this approximation also captures the 
entanglement dynamics, besides some small oscillations 
stemming from neglected rapidly oscillating terms. 



IV. ENTANGLEMENT DYNAMICS FOR 
RESONANT DRIVING 

In order to investigate the entanglement dynamics, we 
consider the exact time-evolution and discuss it within 
RWA. In doing so, we find the value of the effective 
anisotropy parameter 7 determines the qualitative be- 
haviour. 



A. Strong anisotropy 

For a three-spin chain in the limit 7 — )■ 00 ( Jq = 0), the 
repeated action of the Hamiltonian ^ on the initial state 
creates the cyclic sequence |000) |llO) + |011) -> |000). 
This implies that the quantum dynamics is a coherent 
oscillation between only these two states. The corre- 
sponding concurrence reads Ci_3 — | sin(7Ji/2-\/2t)|. In 
particular, at certain times, spins 1 and 3 are fully en- 
tangled, Ci_3 — 1. The exact time-evolution (not shown) 
agrees very well with the RWA prediction. The three- 
spin case also reveals the difference between an open and 
a closed chain: For the closed chain, which is translation- 
invariant, the fully entangled state |110) -I- 1011) would be 



replaced by |110) -I- 1011) -f |101) which has lower bipartite 
concurrence. This emphasizes that lack of translation in- 
variance supports the entanglement creation between the 
ends of the chain. 

For longer chains, the situation becomes more involved, 
but still can be understood qualitatively. Because the 
Hamiltonian conserves parity and the initial state has 
zero excitations, the system will remain at all times in a 
subspace of states having an even number of excitations. 
This together with the fact that the chain is open can 
be used to argue why at resonance there is such a huge 
amount of entanglement. Further, this argument also 
leads to the conclusion that at the point of maximum 
entanglement the reduced state of spins in the ends of 
the chain is (|00) + \n))/V2 for even chains and (|01) -I- 
|10))/\/2 for odd chains, plus a mixed state contribution 
which is smallest the highest the concurrence. A more 
detailed argumentation can be read in the Appendix. 

The resulting entanglement dynamics is shown in 
Fig. [3] We find that the concurrence begins to grow 
after a given time, and reaches a maximum value at 
a time tarrivai- Thereafter, it decays. However, two 
ways of maintaining the achieved concurrence come to 
mind: One can either simply switch off the driving, i.e. 
J{t > iarrivai) ~ Or shift the driving frequency to an 
off-resonant value ujd ^ 2B. The dotted lines in Fig. [3] 
show that both strategies freeze the entanglement as de- 
sired. This certainly requires knowledge of tarrivai which 
behaves very regularly and can be well estimated, as we 
demonstrate below. Moreover, switching off the driving 
parameters has to be much faster than the typical time 
scale of the system, as we assume throughout this work. 
The same appUes also to the onset of the driving. 

B. Moderate anisotropy 

For finite anisotropy 7 ( Jq 7^ 0), the dynamics becomes 
rather complex, see Fig. [3] The concurrence assumes 
several local maxima until the highest one is reached. 
Moreover, we find that the concurrence maximum be- 
came lower. This is due to the presence of swapping 
terms, which spoil the argumentation of the appendix. 
These terms basically will mix the subsets {|00),|11)} 
and {|01),|10)} and thus reduce the maximum achiev- 
able amount of concurrence. 

Thus, we can conclude that the anisotropic limit 7 — > 
00 (Jo = 0) is the optimal working point and, henceforth, 
restrict our discussion to this limit. 



V. SCALABILITY AND ARRIVAL TIME 

Our next goal is to find the arrival time tarrivai and the 
corresponding concurrence maximum as a function of the 
chain length N. Direct integration of the time-dependent 
Schrodinger equation was only possible for up to 12 spins. 
The solution for 8 spins, however, already demonstrates 
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FIG. 4: (color online) Length dependence of the first maxi- 
mum of the concurrence between ends, Ci,]v(t), and the corre- 
sponding fully entangled fraction / for the RWA Hamiltonian 
Q with Jo — 0. The inset shows the length dependence of the 
arrival time at which the entanglement assumes its maximum. 



that for (jjd — 2B, the RWA Hamiltonian ^ captures 
the global behavior very well, see Fig. [3] Therefore, we 
can make even further progress by mapping the RWA 
Hamiltonian to a model for which an exact solution is 
known. For the unitary transformation S = Y[i=i 3 '^f 
which flips the spins with odd site number, we find the 
duality relation 



(6) 



This means that the Hamiltonian for the infinitely 
anisotropic case can be cast as a scaled isotropic Hamilto- 
nian, while our initial state is mapped to the Neel state: 



11010. 



5|0000 . . .). The Hamiltonian H^=q can be 



diagonalized after a Jordan- Wigner transformation [22'. 

Figure [4] shows that the maximum entanglement 
achieved decreases with the chain length rather slowly. 
For very short chains we find almost perfect entangle- 
ment, as predicted above within RWA, while for length 

= 25, the concurrence still possesses the appreciable 
value Ci_25 ~ 0.5. A typical figure of merit in communi- 
cation protocols is the "fully entangled fraction" , defined 
as / = max(e|p|e), where {|e)} is the set of all maximally 
entangled states [53] . Quantum communication protocols 
are superior to their classical counterparts whenever this 
fraction is higher than 2/3. In figure [s] it is seen that 
this magnitude greatly surpasses the classical efficiency 
for rather long chains. 

Note that a chain of length N = 7 represents a partic- 
ular case in which the concurrence equals the fully entan- 
gled fraction. We cannot provide an intuitive explanation 
for this anomalous behavior. 

Already above, we mentioned the importance of know- 
ing the time ^arrival at which the concurrence assumes 
its maximum. In the strongly anisotropic limit 7 — > cx3 
( Jo = 0), we can provide a good estimate for the arrival 
time with the following reasoning: A typical local excita- 
tion will be transported with group velocity Vk = de^/dk, 



where for Jq = 0, the eigenenergies tk = (7J1/2) cos(fc) 
are determined by the wave number k = TTm/{N -f 1), 
m — I, . . . , N and form a band. Since the initial state 
1 0000 . . .) is located in the center of the band, the rele- 
vant wave number is k ~ tt/2. Thus, the time scale for 
traversing the chain is 



t* 



N 



Vn/2 

The inset of Fig. |3] shows that 

^arrival 



2N 
7J1 ' 



1.7 

7J1 



t* 



1.7 + N 

7J1 



(7) 



(8) 



i.e., it grows linearly with the chain length. The factor 
1 /2 on t* reflects the fact that counter-propagating exci- 
tations will meet already in the middle of the chain thus 
establishing distant entanglement. 

Recently, Wichterich and Bose [24 computed the fully 
entangled fraction in spin chains with isotropic nearest- 
neighbor interaction. Starting from the mixed Neel state 
(1/2)10101 . . .)(0101 . . . I + (1/2)|1010 . . .)(1010 ... I, they 
found that switching on a constant interaction entangles 
the spins located at the end sites. The unitary trans- 
formation (|6| maps this model to the limit 7 — > 00 of 
the RWA Hamiltonian Moreover, our discussion of 
the entanglement dynamics within RWA vividly explains 
why in their case the isotropic model permits the creation 
of a remarkably high entanglement. 



VI. IMPLEMENTATION WITH OPTICAL 
LATTICES 

The realization of a XY chain with anisotropy 7=1 
has been proposed for experiments with cold atoms in a 
one-dimensional optical lattice that in transverse direc- 
tion forms a bistable potential [12]. The ground-state 
doublet of each double-well forms the "spin" degree of 
freedom. Then our initial state ^ corresponds to a 
Mott-insulator state, which has already been realized ex- 
perimentally [25J. There the tunnel barriers in longi- 
tudinal direction can be up to ~ 22E^, where the re- 
coil energy typically lies in the kHz regime. This is 
more than sufficient for suppressing longitudinal tunnel- 
ing, such that each double well remains occupied with 
a single atom, while the spin-spin interaction is given 
by a Bose-Hubbard repulsion term. The repulsion term 
is caused by an overlap of Wannier functions describing 
neighbouring atoms. The amount of overlap is given by 
the barrier height, which can be controlled and modu- 
lated via the laser intensity, yielding a time dependent 
J(t). Thus a high barrier effectively yields no overlap 
and hence no spin-spin interaction ( J = 0), whereas a low 
barrier can yield values J ^ 0.1 kHz [26 . The Zeeman 
field B corresponds to the tunnel splitting of the double- 
well potential and is of the order O.IE^ Though, it 
can be manipulated by changing the depth of the dou- 
ble well potential, so that B is greater than, but of the 
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order of, J{t). This implies that the switching times of 
the Zeeman fields have to be considerably smaller than 
1 ms. Coherence times for atoms in such optical lattices 
can be much larger than the system time scale and, thus, 
decoherence should not play a major role. Moreover, 
the initial state |000...0) can be imposed by tayloring the 
field B to be much higher than thermal excitation energy 
ksT due to the environment. Finally, the "spin state" in 
the transverse double well can be probed by fluorescence 
measurement of the atoms. 



VII. CONCLUSIONS 



We have shown that proper ac driving can induce al- 
most perfect entanglement between the first and the last 
spin of an anisotropic XY chain. As a most striking fea- 
ture, we found that the driven chain bears the poten- 
tial for a considerably larger entanglement than the for- 
merly studied static chains. We identified a resonance 
condition which leads to maximal entanglement and also 
provide a reliable estimate for the time after which this 
entanglement is reached. The latter is crucial for freez- 
ing the entanglement once it is created. Our analysis 
within a rotating-wave approximation contributed to a 
qualitative and quantitative understanding of how the 
entanglement is built up: pairwise flipping of neighboring 
spins of an open chain favors correlations between the end 
spins. Moreover, we found that the maximum entangle- 
ment decreases only weakly with the chain length, while 
the entanglement is built up during a time that is lin- 
early length dependent. Thus our protocol demonstrates 
good scalability which is a major requirement for the im- 
plementation of quantum communication protocols. A 
natural application of our scheme is quantum communi- 
cation via state teleportation. This is possible because 
the fully entangled fraction between the first and the last 
spin is sufficiently large, such that a spin singlet can be 
purified |23| . Let us finally emphasize that our protocol 
can be implemented with three different experimental se- 
tups, namely an anisotropic chain with sinusoidal driv- 
ing, an infinitely anistropic chain with a sudden switch 
and an isotropic chain with initial Neel state. This pro- 
vides a broad choice for its application. 
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APPENDIX A: ENTANGLEMENT IN 
STRONGLY ANISOTROPIC CHAINS 

Due to the parity preserving character of the Hamil- 
tonian and the fact that our initial state |00...0) 
has a definite parity, the reduced density matrix of 
the spins at the ends of the chain pn is of the 
form pi|00)(00| +P2(|01)(01| + |10)(10|) -f P3|ll)(ll| + 
(a|00)(ll| + /3|01)(10| + H.c), i.e., in the basis 
{|00),|01),|10),|11)}, it reads 



pi a 

P2 P 

P* P2 

a* P3 



and the corresponding concurrence is 
C = 2max(0, |q;| — P2, | 



(Al) 



(A2) 



Because the chain is open and the Hamiltonian flips spins 
pairwise at adjacent sites, we find (3 = for even chains 
and a = for odd chains. 

For even chains, the proof of this statement is as fol- 
lows: The term |01)(10| stems from evaluating the trace 
over density operators of the form |0[a:]l)(l[j/]0|, where 
the blocks [x] and [y] represent the rest of the chain. Ob- 
viously, only terms with [x] = [y] yield a non-vanishing 
contribution. We demonstrate by reductio ad absurdum 
that it is impossible to fulfill this condition: Let us as- 
sume that states |0[a:]l)(l[a;]0| can occur. By applying 
the pairwise flipping Hamiltonian we have [x\ = [y]l for 
the ket and [x] = 1 [z] for the bra, where the blocks [y] and 
[z\ are yet one spin shorter and, thus, consist of an odd 
number of spins. Hence [x\ = l[a::']l, such that the op- 
erator becomes |01[a;']ll)(ll[a;']10|. Again, by the same 
reasoning we find the requirement [x'\ = 1[?/'] for the 
ket and \x'] = [z']l for the bra. Therefore \x'] = l[a;"]l 
and, thus, |011[a;"]lll)(ll[x']10|. Repeating this proce- 
dure, we end up with a collection of ever smaller blocks 
[a;], [a;'], [x"], [x^"^] all of which possessing an even 
number of spins. Eventually, we remain with the op- 
erator |011...1[a;(")]l...l)(ll...l[a;(")]1...10|. From the ket 
we find the condition that [a;^"-*] = [10], while from the 
bra follows [x^"^] = [01] in order to have a total even 
number of Is and thus conserve parity. Thus we can con- 
clude that initial hypothesis [x\ = [?/] must be wrong. 
This proves that for even chains /? = 0, and so the con- 
currence reduces to C = 2max(0, \ a\ — p^)- This line of 
reasoning can be adapted to the case of odd chains, for 
which one obtains a = 0. 

Yet, in order to obtain a high concurrence, we need 
|ck| 3> P2, as we find in our numerical studies. The trace 
condition for density matrices yields pi -I- p3 = 1 — 2p2, 
while positivity requires a < y^piPa- Note that for pure 
states, a = ^p\Pz ■ Thus, maximizing a necessarily re- 
quires Pi be small, that is, if at any instance of time, a 
starts to increase, as it happens when p2 becomes smaller, 
the concurrence increases as well. Clearly this can occur 
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only at certain times, which is why we see entanglement shared between spins 1 and is then ai|00) +a2|ll) for 
peaks. even chains and ai|01) + a2|10) for odd chains, where we 

At times of maximum concurrence, the resulting state have ignored a small mixed state contribution. 
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